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We analyze the microscopic dynamics of vortex motion through channels that form river-like fractal 
networks in a variety of superconducting samples, and relate it to macroscopic measurable quantities 
such as the power spectrum. As a function of pinning strength, we calculate the fractal dimension, 
tortuosity, and the corresponding voltage noise spectrum. Above a certain pinning strength, a 
remarkable universal drop in both tortuosity and noise power occurs when the vortex motion changes 
from shifting braiding channels to unbraided channels. We compare our results with experiments. 

PACS numbers: 74.60.Ge, 62.20.Fe 



Introduction. — The complex dynamics of a moving su- 
perconducting vortex lattice interacting with material de- 
fects has attracted considerable experimental and theo- 
retical attention, with the observation of intricate chan- 
nels of vortex motion both in simulations , beginning 
with the seminal work in Ref. and through Lorentz 
microscopy Q . Similar channel structures have been ob- 
served in a large variety of systems, including fluid flow 
in a disordered landscape, Josephson junctions, Wigner 
crystals, magnetic bubbles, and stress networks in gran- 
ular systems. These channels resemble the fractal basins 
created by natural rivers |J and other fractal network 
systems (e.g., percolation). 

A quantitative microscopic understanding of the char- 
acteristics of the channels and their effect on macroscopic 
measurements is particularly important for supercon- 
ducting systems, in which the disorder can be controlled. 
Different strengths of disorder produce very different flow 
patterns, ranging from elastic flow to plastic flow [Q, 
characterized by vortices that either remain pinned or 
move intermittently through a vortex river. Since these 
flows can be inferred experimentally via the voltage noise 
produced by the moving vortices |f^-|Io[] , a deeper under- 
standing of the relationship between the noise charac- 
teristics and the properties of the vortex channels would 
lend insight into the experimental systems. 

The complex vortex channel network |]ll"| observed in 
some regimes is difficult to treat analytically. Until now, 
the channel structure has been studied qualitatively in 
simulations , and only transitions caused by changes 
in driving force have been considered (e.g., Refs. 
Transitions in driven vortex lattices caused by differ- 
ent disorder strength (e.g., Ref. Q) are more difficult 
to study since a separate simulation is required for each 
disorder strength. We use a large-scale parallel simula- 
tion to probe 21 samples spanning an order of magnitude 
of pinning strengths, and present a systematic study of 
the transition from one plastic flow phase to another. We 
quantify the fractal nature and the tortuosity of the vor- 
tex channels in the plastic flow for the first time, and 
show how both evolve with disorder strength. We ob- 



serve remarkable correlations between microscopic quan- 
tities such as the tortuosity and macroscopic measures 
such as voltage noise power, corresponding to changes in 
the microscopic nature of the channel flow. 

Simulation. — We model a transverse two-dimensional 
slice (in the x-y plane) of a T — zero-field-cooled super- 
conducting infinite slab containing rigid vortices that are 
parallel to the sample edge (H = Hz). Vortices nucleate 
along one edge of the sample at regular time intervals, 
enter the superconducting slab under the force of their 
mutual repulsion, pass through a pinned region 36A x 36A 
in size (where A is the penetration depth) where a flux 
gradient naturally forms Jl2| , and are removed at the 
other sample edge. Up to fOOO vortices are simultane- 
ously inside the sample, which is periodic only in the y 
direction transverse to the gradient. 

The vortex- vortex repulsion is correctly represented by 
a modified Bessel function, K^r/X). The vortices also 
interact with 1943 non-overlapping attractive parabolic 
wells of radius £ p = 0.15A, representing a density of pin- 
ning sites n p — 1.0/ A 2 . The maximum pinning force, 
fp, of wells in a given sample has a Gaussian distribu- 
tion. We consider 21 samples (a much larger number of 
parameters than in typical simulations |J]) with mean 
values of f p ranging from f p = 0.06/ to f p = l.O/o, 
where f = <&o/87r 2 A 3 . 

The overdamped equation of vortex motion is f; = 
f™ + i" p = r/Vi , where the total force on vortex i 
(due to other vortices f™, and pinning sites f" p ) is given 
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Tik. Here, is the Heav- 



iside step function, (vj) is the location (velocity) of 
the ith vortex, rl is the location of the fcth pinning 
site, £ p is the pinning site radius, N p (N v ) is the num- 
ber of pinning sites (vortices), fy = (r, — ij)/|rj — r,-|, 
= (r, — rj^ )/|r$ — i& |, and we take 77 = 1. We mea- 
sure all forces in units of fo — <&q / 87r 2 A 3 and lengths in 
units of the penetration depth A. We run a highly op- 
timized parallel code on IBM SP parallel computers to 
carefully characterize a wide range of parameters, and we 
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equilibrate each sample for at least 10 6 MD steps before 
taking high resolution data. Due to the open boundary 
conditions in the x direction, the ratio of vortices to pins, 
n v /n p , is not directly controlled. Instead, it emerges nat- 
urally as the system equilibrates. Further simulation de- 
tails appear in Ref. Q. 

Channel network. — We divide the sample into a 300 x 
300 grid to identify the vortex channels. When a vortex 
enters a grid element, the counter associated with that 
grid element is incremented, defining a "channel tran- 
sit" field. All grid elements that are visited at least once 
are considered part of the channel network. We calculate 
the average rate T av at which vortices move through each 
grid site, and construct a distribution of T av over all the 
grid sites to indicate how frequently individual channels 
were traversed. Figure shows channels and distribu- 
tions P(r av ) from four samples with different pinning 
strengths f p . For weak pinning, f p < 0.2/ [Fig. 0(a)], 
the channels cover the entire sample area relatively uni- 
formly. Many grid sites are visited by vortices at a low 
rate, so P(r av ) peaks at small r av . As f p increases, is- 
lands of pinned vortices (shown in white in Fig. 0) form 
and grow. At the same time, favored channels for vor- 
tex flow appear. Grid elements inside channel sites are 
frequently traversed by vortices, so P(r av ) extends to 
higher rates r av . At higher pinning forces when there are 
only a small number of channels, the grid elements inside 
channels produce a distinguishable increase in P(r av ) for 
r a t) ^ 5 x 10~ 4 . The average separation between channels 
d pcrp increases roughly linearly with increasing pinning 
strength f p as the ratio n v /n p increases [Fig. 0(b-d) and 
inset of Fig. |j, until for f p ~ O.66/0, the channels are 
no longer connected in the transverse y-direction, as is 
clearly visible in Fig. 0(d). This breakup represents a 
transition in the nature of the plastic flow, as we shall 
show below. 

Fractal dimension. — To quantify the effect of pinning 
strength f p on the fractal dimension Df of the vor- 
tex channel network, we use a box-counting algorithm 
|p3| to find Df, and plot the results in Fig. 0. Here, 
Df = — lim e _>o log N(e) / log e, where N(e) is the number 
of boxes of side e required to cover all grid sites belonging 
to the channels. The dimension P/ « 2 for low pinning 
strengths, f p < 0.2/o, when vortices are flowing through- 
out the entire sample [Fig. 0(a)]. As f p is increased and 
the channel structure becomes more sparse [Fig. 0(b-c)], 
the fractal dimension Df decreases. 

Our samples with strong disorder have fractal dimen- 
sions close to those predicted recently for channel net- 
works in systems where elastic interactions are unimpor- 
tant [^4 16 1, For example, our sample with f p = 0.75/o 
has a fractal dimension of Df = 1.37, close to the mean- 
field prediction of 4/3 found in Ref. [[uj and the value 
1.38 observed in Ref. [^5). At the strong pinning case of 
f p = 0.83/o, where there are a few isolated channels in 
the sample, we find Df — 1.27. This is in reasonable 
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FIG. 1. Right panels: top view of samples with four dif- 
ferent pinning strengths: (a) f p = O.I8/0, (b) f p = 0.36/o, 
( c ) fp — 0.57/o, (d) f p = 0.66/q. The channels most of- 
ten followed by vortices are darker. The number of channels 
decreases as the pinning strength f p increases. Left panels: 
distribution of the average rate T av at which vortices pass 
through individual grid points in the river. For strong pin- 
ning, the few remaining channels are frequently traveled [see 
brace in (d)]. 



agreement with simulations of non-interacting particles 
@, D f = 1.21, and with other theories Q, Df = 1.22. 
At the very highest pinning strength, f p = 0.9 /q, only 
a single river passes through the sample, giving an ex- 
tremely low fractal dimension Df — 1.15, in agreement 
with the fractal dimension of the main channel of physical 
rivers §,!>/ = 1.14—1.20. 

The fractal dimension gives a static picture of the vor- 
tex channels. We probe the dynamics of the channels by 
considering the fraction N a /N v of vortices that move a 
distance greater than the pinning diameter 2£ p . We find 
that changes in N a /N v closely follow changes in the 
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FIG. 2. Fractal dimension Df (filled diamonds; solid line 
is guide for eye) versus pinning strength f p , and the fraction 
N a /N v (open triangles; dotted line is guide for eye) of vortices 
moving inside the sample versus f p . Both show a slope change 
near f p m 0.6/o- Lower left inset: Average distance between 
channels d pcrp versus f v . Upper right inset: Variation with 
time of the tortuosity, a T , versus f p . Note that ov drastically 
increases near the region where the power So [Fig. 3] peaks. 

fractal dimension Df. At low pinning strengths, f p < 
0.2/ , all of the vortices move, as seen in Figs. 0(a) and||, 
indicating mostly elastic motion. The moving fraction 
N a /N v decreases with increasing pinning strength as the 
motion becomes plastic and some vortices remain perma- 
nently trapped in pinning sites. Fits to N a (f p ) [dashed 
line in Fig. ^] and Df(f p ) [solid line] indicate that there 
is a small but noticeable change of slope in both quanti- 
ties near f p /fo ~ 0.6/o. As we shall see later, this occurs 
when the vortex channels change behavior from braiding 
(f p /fo < 0.6) to non-braiding (f p /f > 0.6). 

Tortuosity and fluctuating braided channels. — To ex- 
amine the motion of individual vortices, we compute the 
tortuosity of the path followed by each vortex. The tor- 
tuosity r measures the amount a path winds [g]: r = x/ L 
where x is the actual distance traveled by the vortex as 
it crosses the sample, and L is the width of the sample. 
Thus, t > 1, and for a straight path, t = 1. 

The plot of the average tortuosity r shown in Fig. || 
reveals a very interesting behavior that is not reflected 
in the fractal dimension Df. For low pinning strengths 
f p < 0.05/o, t ~ 1.1, indicating that the vortex paths 
wind very little. The tortuosity increases with pinning 
force as the pins become more effective and cause individ- 
ual vortex paths to wind around islands of pinned flux, as 
in Fig. |l](b). A vortex can follow a very tortuous trajec- 
tory by crossing between what would have been distinct 
paths at lower pinning strengths. The heavy crossing or 
braiding of channels leads to a peak value of t « 1.5 for 
f p w 0.5/o- As seen in the inset of Fig. |[ the variation 
in time of the tortuosity, o~ T , also peaks near f p ss 0.5/o- 
For f p < 0.5/o, the vortices follow a network of heavily 
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FIG. 3. Tortuosity r (filled diamonds) of vortex paths and 
power So in the second octave of the spectra (open circles), 
both versus pinning force f p . The same results hold for the 
third and fourth octaves. A peak in the noise power was 
also observed in spectra obtained from the voltage trans- 
verse to the driving direction. Left inset: Spectra S(v) for 
f p = 0.18/o. Right inset: Spectra S(v) for f p = 0.66/ . The 
narrow band peaks in both spectra are caused by the regular 
rate at which vortices are added to the sample, while the peak 
at v ~ 10~ 3 in the right inset is due to the typical time of 
flight for motion through the isolated channels. 

braided flow channels. A remarkable drop in r between 
f p ~ 0.5/o and f p ~ 0.7 /o to a saturated low value of 
r ~ 1.15 is clearly visible in Fig. ^. It is important 
to emphasize that the drop coincides with a transition 
from channels that are braided across the entire length 
of the sample at intermediate pinning strengths, 0.2/o < 
f p < 0.5/o [Fig. |l](a-b)], to isolated, unbraided channels 
at higher pinning strengths, f p > 0.6/o [Fig. |l|(d)], that 
are too far apart to significantly interact. This change is 
not merely a finite size effect, since it is occurring over 
length scales significantly smaller than the sample width. 

The crossover from increasing to decreasing tortuos- 
ity results from the combined effects of simultaneously 
increasing the vortex-pin interactions and the flux gra- 
dient. Vortex-pin interactions are less important at low 
pinning forces, and the vortices follow relatively straight 
paths. As the pin strength increases, some vortices be- 
come trapped, the vortex paths begin to wind, and the 
tortuosity increases. The flux gradient is also increasing, 
however, and at the crossover point, the flux gradient be- 
gins to dominate. The vortices then flow directly down 
the steeper gradient, decreasing the tortuosity. We have 
observed the transition in samples of different x direc- 
tion lengths: 18A x 36A, 36A x 36A, and 48A x 36A. The 
pinning force at which the transition occurs shifts down- 
wards slightly as sample length increases. In very long 
samples it is thus possible for the channel flow phase to 
dominate for a fixed current and to be detected by local 
Hall probes. A more detailed account of the effects of 
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sample size will appear elsewhere. 

Voltage Noise. — We next link the transitions in the 
vortex channel structure with experimentally accessible 
voltage noise signals. We sum the forces in the in- 
direction along a strip of the sample 5A in width to obtain 
our voltage signal. We find the spectrum of the result- 
ing signal, V x = J2i £c*^ > f° r eacn sample configuration, 
and indicate the spectral power by plotting the integrated 
noise power in one frequency octave So = dvS{v) ver- 
sus pinning strength f p in Fig. |^. Here, = 1.2 x 10~ 4 
and v\ — 2.4 x 10~ 4 . Remarkably, the form of Sq closely 
follows the tortuosity r. This is because both measures 
are sensitive to the number of metastable states accessi- 
ble to the system. When r is high, the vortices wander 
significantly in the transverse direction, sampling many 
metastable states and thus producing a high noise power. 
The overall drop in noise power for f p > 0.5/o occurs as 
the amount of braiding in the channels decreases. We 
can compare our results to experiments J9|,[l0| , in which a 
peak in the noise power occurs near the depinning transi- 
tion when plastic flow occurs and high r is expected. At 
higher currents, the pinning effectively becomes weaker, 
the vortices flow more elastically, r is expected to be 
lower, and a drop in noise power is observed. This agrees 
with the results shown Fig. ||: for the most plastic flow, 
the highest r occurs and the noise power is highest, while 
for weaker pinning, the vortices flow in straighter paths, r 
is lower, and the noise power drops. Noise measurements 
can thus be used to probe the tortuosity r. 

The shape of the noise curve changes significantly at 
f p ~ 0.5/o [insets of Fig. |). For f p < 0.5/o, the spec- 
trum for v > 1 -3 is of the form S(v) ~ v~ a [left inset of 
Fig. |j, where a decreases for higher pinning strengths. 
For f p < 0.1/o, a « 2, while for 0.2/ < f p < 0.5/ , 
a Fts 1.7, consistent with the experimental measurements 
of D'Anna et al. || and Mar ley et al. ||, respectively. 
During and after the drop in So, for f p > 0.5/o, the spec- 
trum S(v) is no longer of a form that can be characterized 
by a unique slope. Instead, the relatively straight, iso- 
lated channels produce a time of flight signal in the spec- 
trum similar to experimentally observed signals 0. No 
unique time of flight signal appears in our simulation at 
lower pinning forces, f p < 0.5/o, since the many braided 
channels present lead to a spread in the tortuosities and 
a spread in the time spent crossing the sample. 

Summary. — Using novel measures of vortex channel 
structures, including the fractal dimension of the chan- 
nel network and the tortuosity of individual vortex paths, 
we have provided strong evidence for a transition be- 
tween two distinct vortex plastic flow phases as a func- 
tion of disorder strength. We have shown that, as dis- 
order increases, the weak-pinning straight paths contin- 
uously evolve to a braided winding-channel pattern that 
is characterized by high noise power and a high tortu- 
osity with large fluctuations in time. A sharp drop in 



the tortuosity and noise power for intermediate pinning 
signals a transition to flow in non-braiding, isolated in- 
dividual channels. The drop in noise power is consistent 
with experimental measurements ||,[l0|. The transition 
is a universal property of the average pinning strength 
and not a sample-dependent phenomenon, since the dis- 
order configuration was different in every one of the 21 
sample used. This crossover among different dynamical 
flow regimes as a function of pinning strength may also 
be important to other slowly driven disordered systems 
such as Wigner crystals, colloids, Josephson junctions, 
and magnetic bubbles. Our observation that the tortu- 
osity and noise power follow each other closely is a novel 
result indicating that a macroscopic noise power mea- 
surement gives direct insight into the microscopic tortu- 
osity, a link that may be useful for systems with driven 
channels. Our predictions can be tested through experi- 
ments such as Lorentz microscopy or noise measurements 
in superconductors, and direct imaging of colloids. 
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